D-brane Gases and Stabilization of Extra Dimensions in Dilaton Gravity 
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We consider a toy cosmological model with a gas of wrapped _Dp-branes in 10- dimensional dilaton 
gravity compactified on a p-dimensional Ricci flat internal manifold. A consistent generalization 
of the low energy effective field equations in the presence of a conserved brane source coupled to 
dilaton is obtained. It is then shown that the compact dimensions are dynamically stabilized in 
string frame as a result of a balance between negative winding and positive momentum pressures. 
Curiously, when p = 6, i.e. when the observed space is three dimensional, the dilaton becomes a 
constant and stabilization in Einstein frame is also realized. 



One of the main problems of string cosmology is to de- 
termine why the extra compact dimensions evolve differ- 
ently from the observed three dimensions and remained 
comparatively very small. In [1], an intuitive mecha- 
nism was proposed to accommodate this difference where 
strings winding extra dimensions fall out of thermal equi- 
librium and stop the cosmological expansion. In [2-5], 
the arguments of [1] are quantified by demonstrating sta- 
bilization in Einstein and dilaton gravities using the en- 
ergy momentum tensor for string winding and momen- 
tum modes on the torus (see also [6-21] for recent work 
on brane gas cosmology). 

As pointed out in [6], strings may not be able to sta- 
bilize extra dimensions when topology is different than 
a torus. Considering, for example, a spherical internal 
space, there exists no stable winding strings in the spec- 
trum. One would then wonder if the higher dimensional 
branes can play a role in these compactifications. In [6], 
it is shown that a gas of p-branes wrapping over a p- 
dimensional compact, Ricci flat, internal manifold can 
stabilize the volume modulus. As for strings, it turns 
out that there is a balance between the winding and the 
vibrational momentum modes. The results of [6] are ob- 
tained in Einstein gravity and the purpose of this letter is 
to generalize the framework to dilaton gravity. Here, it is 
natural to consider a toy model with a gas of winding D- 
branes. Since dilaton couples to Z?-brane world-volume, 
it should be activated in this scenario. 

Our first aim is to generalize the equations of motion 
in the presence of a conserved source coupled to dilaton. 
In string frame, the low energy effective action (we set 
H^_,up = 0) is given by 



field equations can be found as 



j S^'x^e-^^ [R + 4(V0)2 + e^^Cr,,] , (1 



where we add a term to include the effects of matter 
governed by the Lagrangian Cm- From this action, the 
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where Tnn is the matter energy momentum tensor 
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and the extra contribution F in (3) arises from the cou- 
pling of dilaton to Cm- 

The conservation formula V^T'*" = plays the role of 
matter field equations. To determine the unknown func- 
tion F we note the well known fact that (3) can be viewed 
as a consequence of (2) and the contracted Bianchi iden- 
tity. Taking the divergence of (2) with V^, one finds: 



(5) 



Although in general F cannot be uniquely fixed by this 
constraint, in a cosmological context one can assume that 



0, 



(6) 



where t is the time coordinate and i denotes a spatial 
direction. Under these conditions, (5) gives 

F = (a - 2) e'"^ Tu g" = - (a - 2) e°-^ p. (7) 

It is easy to see that (7) corresponds to the choice 

C,n = -p, (8) 

which is precisely the Lagrangian for hydrodynamical 
matter (see e.g. section 10.2 of [22]). 

Let us take a (1 + to + p)-dimensional space-time with 
the metric 



-dt^ 



dx'dx' + R^ dS^, 



(9) 
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and consider a p-brane wrapping over Ep. Here i, j = 
l,..,m, (a, 6 = Sp is a p-dimensional compact 

manifold and R is a constant scale factor. The brane 
dynamics is governed by the usual action 



(10) 



where Tp is the tension, 7 is the induced metric and the 
local brane coordinates are identified as a = {t, Sp). It is 
easy to see that the wrapped brane ofTcrs a stable classi- 
cal configuration which solves the embedding equations 
derived from (10). Although it is not known how to quan- 
tize (10) exactly, for our purposes it is enough to employ 
an approximate scheme and simply consider small fluctu- 
ations around the classical background. It is clear that in 
this theory there are both winding and vibrational mo- 
mentum modes in the spectrum. The energy of a winding 
mode is given by 



E = {nTpQ.p)BF, 



(11) 



where n is the winding number and fip is the volume of 
Sp. The spectrum of small vibrations can be found by ex- 
panding the action (10) to the second order in transverse 
fiuctuation fields which yields free modes propagating in 
Sp. The energy spectrum is given by [6] 



i?' 



(12) 



where — A| are the discrete eigenvalues of the Laplacian 
on Sp. For D-branes, there are additional U{1) gauge 
fields living on the world- volume. At the linearized level, 
the transverse excitations obey free field equations and 
the energy spectrum is identical to (12). 

The average pressure of a closed system can be foimd 
from the energy as a function of volume; EiV). Recalling 
that in our case the volume modulus is fixed by BF one 
gets 



P = 



-(n Tp rip) : winding, 
Afe/(pi?) 



momentum. 



(13) 



Note that the winding and the momentum modes apply 
negative and positive pressures, respectively. For mo- 
mentum modes, 1/R dependence of the energy and the 
pressure can also be inferred from the uncertainty prin- 
ciple and the fact that the size of a compact direction is 
approximately equal to R. 

For a non-interacting gas of such excitations, it is 
enough to sum over the modes with constant number 
densities to find the total energy and pressure. However, 
it is more natural to assume a state of thermal equilib- 
rium^ and consequently allow transitions between differ- 
ent energy levels. In that case, the mean number of a 



mode in the gas is proportional to the Boltzmann factor 
e~^^ . At very low temperatures the gas mainly consists 
of the lowest energy modes since the relative densities of 
other excitations are exponentially suppressed. In this 
case, one can use the expressions (11)-(13) with fixed 
quantum numbers (corresponding to lowest lying states) 
to obtain the average energy and pressure of the gas. 
In a cosmological setting, one should replace (9) with 

ds^ = -dt^ + e'^^dx^dx' + e^^ d^l, (14) 

where the metric functions B and C depend only on the 
proper time t. In this work we assume that dT,^ is Ricci 
flat. The scale factors for the observed and the internal 
spaces are defined by 



Rob — e 



Ri'. 



(15) 



To find the energy momentum tensor from (11)-(13), the 
constant R in these formulas should be replaced with the 
time dependent function e'~' and one should also divide 
by the total spatial volume 6"*^+^*^ to find the corre- 
sponding densities. This gives 

iB , ^-mB-{p+l)C 



— e 
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(16) 
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where and Tm are constants. It is easy to check that 

T^^ is conserved; V^.T'"" = 0. As noted in [6], (16) 
implies an equation of state = Wj p, where for winding 
modes 



— 1 : brane direction, 
: transverse direction. 



and for momentum modes 



1 /p : brane direction. 



(17) 



(18) 



: transverse direction, 

which is equivalent to radiation confined in the compact 
space. 

Up to now our considerations are general and the for- 
mulas are valid for diff'erent brane types. To focus on 
D-branes we should fix the value of a. Since Z)-branes 
contribute first at disk order in string perturbation the- 
ory, their tension is inversely proportional to string cou- 
pling Qs = e'^. Therefore, one should set^ 



a= 1 



(19) 



^ As discussed in [4], gravitons and photons offer an ideal can- 



didate for a thermal bath which does not affect the dynamics 
of extra dimensions. In [6] we also observe that brane winding 
and momentum modes dominate the cosmology at late times in 
the presence of matter which has equation of state p = ujp with 
> 0. Therefore, the contribution of the thermal bath can be 
ignored in the following. 
^ For fundamental strings the tension does not depend on and 
one should set a = 2. In that case the field equations (2)-(3) 
agree with that of [2]. 
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to have the correct Qs dependence in (1). Recall that we 

have also 



m + p = 9 



(20) 



in string theory. 

Using (16) in (2), (3) and (7) one can get the following 
equations (we take some linear combinations to simplify 
the expressions): 



3-p 



Fi +2F2, 



(21) 



B= -kB+^F^ + ]^F2, 

C= -kC-\F^+^F2, 
2 2p 

= mB^ + pC'^ + 2Fi + 2F2, (22) 



where 



k = mB + pC - 2^, 

Fi = T^e'^-'"^, 

- T p'l>-mB-{p+l)C 



(23) 



and dot denotes differentiation with respect to t. Eq. 
(22) can be viewed as a constraint on initial data and 
when it is obeyed at some to it will be fulfilled during 
subsequent evolution. 

We could not obtain the most general solution of these 
equations. Using, however, the following ansatz 



B = bi ln{t), 
C = Co, 



(24) 



we find that (21) gives 
6.= " 



2 (m - 3) 

T^' '^1 ~ To"' 

m + 3 m + 3 

A{p + 2)p 



T„(p+l)(m + 3)2^ 
(p+l)Co = {P + '2)Tm 



(25) 
(26) 
(27) 



and (22) is identically satisfied. 

In summary, we have the following power-law solution: 



ds^ = -dt" + {tfl^^+^) dx'dx' + e^^" dEl 



^<t> = g<^o (■^^2(m-3)/(m+3) 



(28) 



It is obvious that in this background the internal dimen- 
sions are stabilized in string frame. Remarkably, when 
the observed space is three dimensional (i.e. when m = 2>) 
the dilaton does not depend on t and the canonical Ein- 
stein frame becomes identical to string frame. Again for 
m = 3, the expansion of the observed space is equivalent 



to the one for pressureless matter in standard cosmology. 
When TO > 3 the dilaton grows and at some point we 
enter the strong coupling regime in which the low energy 
field equations can no longer be trusted. 

Apparently, (28) corresponds to special initial condi- 
tions and to be able to talk about dynamical stabilization 
we should consider more general initial data. To verify 
stabilization we did several numerical integrations with 
arbitrary initial conditions. Before presenting our find- 
ings let us discuss the following issue. The constraint 
(22) divides the solution space into two disjoint pieces 
according to the sign of k defined in (23), since by (22) 
k cannot vanish and thus its sign cannot alter in time. 
The sign is strongly correlated with the arrow of time 
in the problem. In all our numerical integrations we ob- 
serve that when fc < a singularity at finite proper time 
is encountered. These solutions can be thought to evolve 
backwards in time and hit the initial "big bang" singu- 
larity in finite duration. On the other hand all numerical 
integrations with > yield smooth solutions and these 
can be thought to evolve forwards in time. 

There is also a physically unexpected behavior for neg- 
ative k. From the second equation in (21) we find that 
for fc < and when initially B > 0, B is always positive^ 
which implies an accelerated expansion for the observed 
space since the scale factor consequently obeys Rob > 0. 
This, however, contradicts with the anticipation that the 
winding and the momentum modes produce a decelera- 
tion by slowing down the expansion since from (17) and 
(18) they do not apply any pressure along the observed 
directions, i.e. they act like pressureless dust. This shows 
that (at least some) solutions with A; < are not physi- 
cally well-behaved. 

Focusing, therefore, on initial data with fc > 0, we ob- 
serve that in all numerical integrations the metric and 
the dilaton asymptotically approach (28). In figures 1, 2 
and 3 we plot two illustrative examples for m = 3, 2, 4, 
respectively. In these runs while B eventually increases 
in time, C performs damped oscillations around the con- 
stant value (27). For to > 3 the dilaton grows, for to < 3 
it decreases and for to = 3 it also oscillates around (26) 
which are damped in time. These results strongly indi- 
cate that (28) is a future asymptotic "fixed point" for 
(21) and (22). 

To support these numerical findings analytically let us 
note that (21) can be viewed to describe a point particle 
in (4), B, C) space. There are three forces acting on the 
particle: two of them are characterized by the functions 
Fi^2 and there is also a frictional force which is deter- 
mined by the time dependent coefficient k. The fact that 
we restrict /c > is crucial in identifying the frictional 



^ It is easy to see that the analogous argument for fe > does not 
work. 
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FIG. 1: (m,p) = (3,6): The graphs of {B,C,(t)) for the fol- 
lowing initial conditions = (0.1, 2.0), B(l) = (-2.0, 0.3), 
(7(1) = (0.1,-0.3), (7(1) = (1.7,0.4) 0(1) = (0.2,0.7), 
^(1) = ( — 1.0,1.0). respectively. We choose T^, = 4 and 
Tm = 3 so that the internal dimensions are stabilized at 
(7 = 0. 



force*. Ignoring friction for a moment, the third equatfon 

in (21) shows that the driving force along C vanishes at 
C = Co, where Cq is given by (27). (7 = (7o is a stable 
equihbrium point since when C > Co the net force is neg- 
ative and when C < Co it is positive. For stabilization it 
is also crucial that the oscillations around Co are damped 
and this is exactly achieved by friction. Therefore, as 
t ^ 00 one expects C ^ Co which proves stabilization. 

On the other hand, using the asymptotic value Co in 
(21) one sees that the sum of the forces generated by 
Fi_2 in (f) direction only vanishes for m = 3 {p = 6). In 
this case, due to the existence of friction the motion in (f) 
should eventually stop and this explains the stabilization 
of dilaton when m = 3. For m 7^ 3, the net force never 
vanishes and thus will be time dependent. Similarly, 
the second equation in (21) shows that the net force along 



* For fc < there is an "anti-frictional" force. Not surprisingly 
this would destabilize the system and create a singularity. The 
existence of friction rather than anti-friction also signals that 
flow of time in the problem is correctly identified. 



B 




FIG. 2: (m,p) = (2,7): The evolution of {B,C,<I)) for the 
initial data 5(1) = (0.7,2.0), B(l) = (0.9,-0.5), C(l) = 
(-0.2,0.2), (7(1) = (1.7,-0.3) 0(1) = (2,-0.7), 0(1) = 
(0.3, —2.1), respectively. We set T^, = 9, Tm = 7 and as a 
result (7o = from (27). 

B direction is never zero and thus B cannot become a 
constant. 

Technically, it is necessary for stabilization that Fi and 
F2 terms in the third equation in (21) have opposite signs. 
The value of a is crucial for this condition to hold. In 
general the equation for C takes the form 

C=-kC-lF, + '-±Sl^F2. (29) 
2 2p 

For a < both terms become positive and stabilization 
cannot be achieved. In the interval < a < 2 the forces 
have opposite signs and for a > {2/p) + 2 the stabilization 
is again ruined. Remarkably for _D-branes a = 1, which 
is in the range that allows stabilization. 

Concluding, in this paper we have presented a dynam- 
ical stabilization mechanism for the volume modulus in 
string theory based on the cosmological impact of a gas 
of £>-branes. The model is obviously incomplete, how- 
ever it reveals a generic dynamical behavior. It is not 
clear if this mechanism can be used in a realistic sce- 
nario to ensure moduli stabilization. One problem is 
that a straightforward generalization of Brandenberger 
and Vafa (BV) mechanism [7] implies all p-branes with 
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FIG. 3: {m,p) = (4,5): The results of numerical integra- 
tions with the initial conditions B{1) — (1.0,0.3), B{1) = 
(-1.0,0.1), C(l) = (0.2,-0.1), C(l) = (0.0,0.2), ,^(1) = 
(1,-0.7), 0(1) = (-3.0,-0.4), respectively. We have also 
T^=7 and = 5. 



p > 2 to annihilate in the early universe. Therefore, 
having a gas of higher dimensional branes at late times 
seems unlikely. On the other hand, as emphasized be- 
fore, it seems that higher dimensional branes are needed 
for moduli stabilization in compactifications on topolog- 
ically non-toroidal spaces like spheres or phenomenolog- 
ically viable Calabi-Yau or G2 manifolds. Therefore, to 
have a realistic scenario involving brane gases it seems 
that we should either modify BV mechanism or improve 
the method based on winding strings to circumvent the 
topological restrictions. 
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